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Abstract. This review paper is devoted to a presentation of rceent progress in 
wave turbulence. I first present the context and state of the art of this field of 
research both experimentally and theoretically. I then focus on the case of wave 
turbulence on the surface of a fluid, and I discuss the main results obtained by 
our group: caracterization of the gravity and capillary wave turbulence regimes, 
the first observation of intermittency in wave turbulence, the occurrence of 
strong fluctuations of injected power in the fluid, the observation of a pure 
capillary wave turbulence in low gravity environment and the observation of 
magnetic wave turbulence on the surface of a ferrofluid. Finally, open questions 
in wave turbulence arc discussed. 

1. Context and state of the art. When strong enough surface waves propagate 
in a medium, their interactions can generate waves of different wavelengths. This 
energy transfer through the different spatial scales can occur on a wide range of 
wavelengths. This stationary out-of-equilibrium state is called the wave turbulence: 
the energy of the system cascades through the scales from a scale where energy 
is injected up to a small scale where energy is dissipated. Wave turbulence thus 
concerns the study of dynamical and statistical properties of an ensemble of in- 
teracting nonlinear waves. The archetype of wave turbulence is the study of the 
random state of ocean surface waves generated by wind or current. But this is a 
very common phenomenon that is present in various situations: surface waves on 
the sea [28, 42, 50, 90], internal waves in the ocean [64], Alfven waves in solar winds 
[82], radar waves in ionosphere, spin waves in solids, Rossby waves in geophysics, 
ion waves [66] and Langmuir waves [49, 96] in plasmas, waves in nonlinear optics 
[59], quantum waves in Bose condensates [9, 24, 30, 80] ... 

Wave turbulence is thus an interdisciplinary subject that involves different com- 
munity: astrophysics, geophysics, mathematics, fluid mechanics and different fields 
of physics. Since the early developments of theoretical tools by applied mathe- 
maticians and physicists, the first ones to be interested on wave turbulence were 
occanographcrs and meteorologists. Their motivations are numerous such as to de- 
velop climatic models, predict the sea state with more precision, or extract some 
energy from ocean waves as a source of alternative energy... At a more fundamental 
level, the goal is to understand the energy transfers between nonlinear interacting 
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waves by means of generic laws whatever the particular medium in which these 
waves propagate. 

Contrarily to a widespread belief, the analogy between wave turbulence and 
hydrodynamic turbulence within a fluid (3D turbulence) or within a film of fluid (2D 
turbulence) is quite limited. Although being also governed by the nonlinear effects 
and studied in a statistical and non-determinist way, wave turbulence is described 
by equations which strongly differs from Navier-Stokes ones of usual turbulence. 

Wave turbulence theory, also known as weak turbulence (WT), is a statistical 
theory describing an ensemble of weakly nonlinear interacting waves. The first ana- 
lytical studies began in the sixties [7, 8, 46]. They were motivated by oceanographic 
measurements of Fourier spectra of the surface wave elevation as a function of the 
wavenumber. The aim was to understand how nonlinear interactions between the 
waves generate transfer of energy among different scales thus explaining the shapes 
of the observed Fourier spectra. In the early seventies, Zakharov and co-workers 
understood that, besides equilibrium spectra, nonlinear ly interacting waves also in- 
volve spectra that correspond to the transfer of a finite energy flux from the large 
scales to small ones [97, 98, 99]. These stationary out-of-equilibrium solutions car- 
rying fluxes of conserved quantities are the analogue of the Kolmogorov spectra of 
hydrodynamic turbulence. However, in the case of interacting waves, they can be 
computed analytically. To wit, let us consider a wave with a wavenumber k and a 
frequency w^. Its energy can be read, 

% = nj: oJk (1) 

where is called the "wave action". WT then leads to the expression of the 
temporal evolution of the spectral density of the wave action by means of a kinetic 
equation such as 

^ = + (2) 

where Cj: is an interaction integral between waves {e.g. analogue of the collision 
integral of the Boltzmann formalism of the kinetic theory of gases) , Z?^ a damping 
term representing the energy dissipation and a forcing term describing the energy 
injection which generates the waves. When the nonlinearity in is weak and the 
number of excited modes is large, a statistical description of wave turbulence makes 
sense. Indeed, by means of ensemble averages, the combination of weak nonlinearity 
and dispersion ensures that the long time asymptotics of the probability distribution 
of the wave field are universal (close to a Gaussian) for a broad band of initial 
distributions. This means, that over sufficiently long times, the system relaxes 
towards Gaussian statistics insensitive to the initial probability distribution that is 
unkown a priori. This is called the asymptotic closure. The stationary solutions 
of the kinetic equations of weak turbulence can be then computed analytically at 
the equilibrium or in an out-of-equilibrium regime using perturbation methods of 
a small parameter quantifying the strength of the nonlinear term relative to the 
linear one. The stationary solutions are found scale invariant for isotropic and 
homogeneous systems by using a set of conformal transformations [99]. At the 
equilibrium {Dj: = = 0), the classical solution is the so-called Rayleigh- Jeans 
spectrum corresponding to the equipartition of energy between modes. An out-of- 
equilibrium solution {Dj: 7^ 0, J^: ^ 0) is tractable when the scale of the energy 
injection (large scale) and the scale of the dissipation (small scale) are assumed 
widely separated to obtain an inertial regime independent of both the injection and 
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the dissipation. Eq. (2) then has conserved quantities (energy flux e = J n^cjkdk 
or "particle" flux). For the finite-energy flux solution, the non zero first order 
expansion of the interaction integral allows to exactly works out the spectral density 
of the wave action nj: or the spectral density of the energy Ej: such as^ 

% ^ e^/^^-^^A:"" with a > f , (3) 

where the exponent of the energy flux depends on the N-wavc interaction process 
which is itself fixed by both the wave dispersion relation and the dominant nonlinear 
interaction^. This out-of-equilibrium solution is called the Kolmogorov-Zakharov 
spectrum (by analogy with the power-law spectrum obtained dimensionally for usual 
turbulence). Consequently, a direct "energy cascade" occurs: due to the interac- 
tions between nonlinear waves, the energy flux injected at large scale is transfers 
towards smaller scales. Similarly, for the finite-particle flux solution, the wave ac- 
tion J\f = J n-j^dk is conserved, an invariant scale solution is found and an inverse 
cascade occurs (from small scales to large scales)''. These spectra have been exactly 
computed in nearly all fields of physics involving wave dynamics [99] . A lot of data 
obtained by remote sensing of the atmosphere or the ocean as well as satellite mea- 
surements in astrophysics, have thus been analyzed using the framework of wave 
turbulence. 

The main assumptions for the derivation of the weak turbulence theory are: 

• A weak nonlinearity is required due to the perturbation methods: the non- 
linear term has to be much smaller than linear one (the fast linear time scale 
is separated from the slow nonlinear time scale). In practice, this condition 
corresponds to waves with enough amplitudes but not too high. 

• The size of the system is assumed infinite. No boundary effect is taken into 
account theoretically. 

• The conserved quantity (energy flux) in wave turbulence is assumed constant 
with no fluctuation during the cascade through the scales. 

• The hamiltonian formulation of weak turbulence assumes no dissipation at 
least in a transparency window in between the scale of injection (usually as- 
sumed at A: = or at the scale of the system), and the scale of dissipation at 
k — > oo. In practice, this means that the forcing of waves has not to be at all 
scales. 

• The system is supposed isotropic and homogeneous spatially, and no coherent 
structures'* are involved. 

• Waves have random phases. This conditions allows to derive a statistical 
equation for the ensemble averaged correlations of wave amplitudes. 

• The single assumption on the initial distribution of the wave field is that its 
cumulants have to decay at infinity. 

At the moment, theoretical works try to model some finite size effects [51, 52, 53, 
70, 88, 100]. The role of coherent structures on the dynamics also begins to take 
into account theoretically [23]. During the last decade, theoretical or numerical 

^Note that, for isotropic systems, a frequency power law spectrum corresponds to the spatial 
one of Eq. (3) by using the dispersion relation. 

^For instance, the first non-zero resonant term for gravity waves is a = 4 wave interaction 
process, whereas for capillary waves N = 3. 

•^This is the analogue, in 2D turbulence, of the enstrophy conservation for the direct cascade, 
and the energy conservation for the indirect cascade. 

*For example, in hydrodynamics, vortices, cusps or wave breakings. 
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works go beyond spectrum prediction. Transient behaviors during the evolution 
towards the stationary regime (or decaying wave turbulence) [55, 56], condensation 
of classical waves (and its analogy with Bose condensation) [9, 24, 30, 58, 80] has 
been performed as well as the introduction of an empirical dissipation to model 
some strongly nonlinear processes (such as wave breakings) [101]. 

In spite of the numerous analytical predictions for more than 40 years, wave 
turbulence is not so much studied experimentally. This is thus the opposite situation 
of the usual turbulence in which numerous experimental results exist, but very few 
analytical results can be carefully derived from master equations. 

Oceanography and atmospheric sciences provide more and more in situ data 
nowadays given by weather balloons, surface buoys, radar and satellite measure- 
ments [72, 92]. However, these large-scale systems depend on numerous parameters 
(namely for the ocean surface waves: wind, oceanic current, fetch...) leading to 
empirical determinations of the energy spectral density [72]. Wind on the ocean 
surface also generates a wave forcing at various scales, and thus without scale sep- 
aration as imposed by the theory. Laboratory experiments are much more relevant 
to accurately tune and control the system parameters. Surprisingly, such laboratory 
experiments are rare, but are fundamental in order to compare them with in situ 
measurements and theoretical predictions of weak turbulence. Only few laboratory 
experiments have been performed since 1996, mostly on wave turbulence on the 
surface of a fluid [14, 15, 47, 48, 62, 93, 94]. These works were focused on the 
measurement of the frequency power spectrum of capillary waves in roughly good 
agreement with the theoretical prediction of weak turbulence [98]. Note that, some 
recent experiments of elastic wave turbulence on the surface of a thin plate [10, 68] 
show that the energy spectrum is in strong disagreement with the predictions of 
weak turbulence [29] . 

In this paper, I will focus on the hydrodynamic case of wave turbulence. In § 2, 
I will first review the previous results on this issue. Then, the main results of our 
experimental works on wave turbulence on the surface of a fluid will be described 
in § 3. Finally, § 4 will be devoted to the conclusion and a discussion on open 
questions in wave turbulence. 



2. Hydrodynamic v^rave turbulence. Waves on the surface of a fluid are in a 
gravity regime when their wavelength A is large with respect to the capillary length 
Ac = 2nlc = 2ti^J^ / pg where g is the acceleration of gravity, h the fluid depth, p 
the density and 7 the surface tension. Otherwise, they are in a capillary regime. 
The crossover between gravity and capillary regimes occurs for a wavelength close 
to the centimeter whatever the working fluid. The dispersion relation of linear 
inviscid surface waves reads uj^{k) = \gk + ["f / p)k'^'\iai\\i{kh). In a deep fluid 
approximation, when A <C 27r/i, one has 

= gk +-k^ (4) 
P 

By balancing both terms of the second hand of Eq. (4), one finds the critical 
wavenumber kc = 1/lc for the transition between both regimes. If waves are weakly 
nonlinear, the weak turbulence theory leads to the expression of the power spectrum 
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of surface wave amplitude such as 




for gravity waves [97, 103] 



for capillary waves [98] 



(5) 



(6) 



These relations can be also computed by dimensional analysis. In the gravity 
regime, parameters involved are the power spectrum of wave amplitude 5^, the 
gravity g, the energy flux e and the frequency u), of dimensions: [S'^] = i^T ; 
[e] = L^T~^ ; [g] = LT~'^ ; [w] = T"-*^. One can conclude by dimensional analysis 
only if the scaling law between Sf^^"''" and e is assumed. As explained above in §1, 
this scaling is known as soon as the number TV of interacting waves is known [i.e, 
the order of the first non-zero term in the interaction integral Cj: of Eq. (2)]. Thus, 
for a N-wave resonant process, the power spectrum scales as £^/(^~^). Gravity 
waves involving a 4- wave interaction process, one has S^"^""^ oc e^/"^. One thus finds 
dimensionally the power spectrum of gravity wave as Eq. (5). Similar dimensional 
analysis can be performed for capillary waves assuming a 3-wave interaction pro- 
cess. Note that this dimensional analysis can be done whatever the wave dispersion 
relation of the system [23]. 

Ocean surface waves are mainly generated by winds and currents at various scales. 
In the literature, in situ measurements of the exponent of the frequency-power law 
of the wave amplitude spectrum varies considerably, even if certain measurements 
roughly leads to a uj~^ scaling [28, 42, 50, 90] thus suggesting a possible agree- 
ment with weak turbulence theory [97, 103] and numerical simulations [31, 73, 78] 
of gravity wave turbulence. However, the wave spectrum depends on numerous 
uncontrolled parameters such as duration of wind blowing, fetch length, stage of 
storm growth and decay, and existence of a swell, leading to fit the spectrum with 
too many parameters [72]. 

When gravity wave turbulence is not forced by wind, mechanisms of energy 
cascade in the inertial regime change, as well as the scaling law of the spectrum 
[54, 60, 73]. Indeed, if does not depend on the energy flux, Phillips has shown 
dimensionally since 1958 that Srj{u!) oc £^g'^uj~^ [75]. This situation corresponds to 
gravity waves of high amplitudes dissipating all the injected power (by wave break- 
ings for instance). More recently, Kuznetsov has taken into account the presence 
of possible sharp wave-crests (cusps) on the surface that are assumed to propagate 
without deformation (i.e. using w cx fc instead of a; = \/gk). If these sharp-crested 
structures occur along ridges (lines) then Srj{uj) oc w"'* [60]. This exponent is similar 
to the one in Eq. (5) computed by weak turbulence theory (where no crested waves 
are involved). In the same way, when these wave slope divergences are assumed to 
be isolated peaks (points) propagating as a; = \/gk, the Phillips spectrum is found 
again [75]. Recently, theoretical and numerical attempts try to take into account 
the quantization of wavenumber k in finite size systems that can strongly affect 
energy transfers [51, 52, 53, 70, 88, 100]. Notably, for gravity waves in a tank of 
characteristic size L, Nazarenko predicts that S'^(w) oc e'^ g^l'^L~^l'^ijj~^ [70]. 

Experimentally, the capillary wave turbulence regime has been observed by means 
of optical methods [14, 15, 47, 48, 62, 93, 94] and with a parametric forcing of the 
waves by shaking the whole container holding the fluid. It has been reported that 
the height of the surface displays a power-law frequency spectrum in w^^^/^ in 
agreement with weak turbulence theory [98] and simulations [79]. However, peaks 
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Figure 1. Photo and schematic view of the experimental setup. 
Wave makers arc on both sides of the capacitive sensor measuring 
the wave amplitude at a given location. Tank size: 20 cm x 20 cm 
filled with water or mercury. 

and its harmonics (related to the parametric forcing) arc observed on the spectrum 
with maximal amplitudes decreasing as a frequency power-law [14, 15, 87, 93, 94]. 
The frequency-spectrum exponent estimated in this way is thus not very accurate. 
A recent study has even found an exponent in disagreement with weak turbulence 
underlying the difficulty to reach a wave turbulence regime with a parametric forcing 
[87]. Finally, note that new experiments of gravity wave turbulence in a large tank 
have been recently carried out [27, 83, 84]. 

3. Overview of our results. We have developed experiments to study wave tur- 
bulence on the surface of a fluid forced by wave makers [11, 34, 35, 36, 38]. This 
investigation is of particular relevance since it allows to study on the same system 
two different regimes, those of gravity and capillary waves, that are characterized 
by different nonlinearities, 4-wave and 3-wave interactions respectively, leading to 
different poweer-law spectra. The crossover between the two regimes is also inter- 
esting, since it presents deviations from the weak turbulence theory"". The role of 
the forcing parameters on the spectra had never been reported, nor the existence of 
a possible intermittency. The energy flux and its fluctuations had never been mea- 
sured in wave turbulence, although it is the conserved quantity in wave turbulence. 
Another motivation is also to understand the statistical properties of energy flux 
necessary to maintain a dissipativc system in an out-of-equilibrium stationary state. 
This is a general issue that is not restricted to the single case of wave turbulence. 

3.1. Spectrum and probability distribution of v^fave amplitudes [34]. As 

shown in Fig. 1, surface waves are generated by two wave makers driven with a 
random excitation (both in frequency and amplitude) by means of electromagnetic 
vibration exciters. This random forcing is in a narrow low-frequency range selected 
by a low-pass filter (typically 0.1 - 5 Hz). The square tank, 20 cm side, is filled 
either with water or mercury up to a height of 2 cm®. The wave amplitudes are 
measured at a given location by means of a capacitive method. This energy injected 
to the system by large wavelengths is transferred towards small structures due to 

^ Close to this transition, the nonhncarity is not weak and the dynamics is dominated by fully 
nonlinear structures [23]. 

®Wavcs of wavelengths A <C 2'iTh ~ 12 cm are thus in a deep-water regime. 
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Figure 2. Typical time recording of the surface waves, r]{t), dur- 
ing 50 s. {t]) = 0. 




Figure 3. Power spectrum of wave amplitude showing the regimes 
of gravity and capillary wave turbulence. The transition ~ 17 Hz 
between both regimes corresponds to Ac = 2ttIc — 1 cm, where Ic 
is the capillary length. Dashed lines have slopes —6.1 and —3.2. 
Random forcing with a narrow bandwidth < 4 Hz. 



nonlinearities. This process is characterized by measuring the frequency spectrum 
and the distributions of wave amplitude fluctuations. 

A typical recording of the surface wave amplitude, ri(t), is displayed in Fig. 2 
as a function of time. The signal is very erratic and the asymmetry observed is a 
signature of the steepness of the waves: high crest waves are more probable than 
deep trough waves. By Fourier transform of such a signal, one can obtain the power 
spectrum density of the wave amplitudes 5*^(0;) = / + T)rj{t)) e~'^'^'^ dr . Waves 
being in nonlinear interactions, a wave turbulence regime is thus expected, that is 
the observation of a scale invariant spectrum. 

As shown in Fig. 3, we have observed, the crossover between the regimes of gravity 
(at large scales) and capillary (at small scales) wave turbulence [34]. Spectrum 
displays two different power laws corresponding to each regime. In the capillary 
regime, the frequency power-law exponent is in agreement with weak turbulence 
theory of Eq. (6) in w^^^/^. Our estimation of the capillary-wave spectrum exponent 
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Figure 4. Slopes of spectra of capillary (open symbols) and grav- 
ity (solid symbols) wave turbulence for different forcing amplitudes, 
cr^, and bandwidths [(o) to 4 Hz, (y) to 5 Hz and (□) to 6 
Hz]. Dashed lines are weak turbulence predictions [see Eqs. (5) & 
(6)]. Inset: Crossover frequency between both regimes as a function 
of the forcing parameters. 




Figure 5. Inset: Probability density functions of the wave ampli- 
tude, T]{t), for different increasing amplitude of the forcing (from 
black to red curves). Main: Same distributions rescaled by its stan- 
dard deviation yijf). Dashed lines: Gaussian fit with zero mean 
and unit standard deviation. 



is much more accurate than in experiments using a parametric and monochromatic 
forcing [14, 15, 47, 48, 62, 93, 94]. The transition between both regimes corresponds 
to a wavenumber of the order of the inverse of the capillary length l^. For usual 
fluids, this corresponds to a critical frequency of fc — \/ g/2lc/'K ~ 17 Hz, that is 
a wavelength of the order of 1 cm. In the gravity regime, the spectrum exponent 
is found to be dependent on the forcing parameters (amplitude and bandwidth of 
the random forcing) as shown in Fig. 4, and consequently in disagreement with the 
weak turbulence prediction of Eq. (5) in The crossover frequency between 

both regimes also depends on the forcing parameters (see inset of Fig. 4). This 
dependence of the exponent of the gravity wave spectrum has been recently found 



LABORATORY EXPERIMENTS ON WAVE TURBULENCE 



9 



again in a much larger tank [27]. Its origin is still an open problem. It could 
stem from finite size effects of the tank [100]. It has been also shown numerically 
that gravity wave spectrum is very sensitive to the condensation of long waves 
background (or "condensate" following analogy with Bose-Einstein condensation in 
condensed matter physics) [9, 24, 30, 58, 80]. 

At low forcing amplitude, the probability density function (PDF) of the wave 
amplitudes is found to be Gaussian (see Fig. 5). At high enough forcing, the PDF 
becomes asymmetric thus showing the non-Gaussian feature of wave turbulence 
[34]. This phenomenon is well known by oceanographers which measure crest-to- 
trough wave distributions which deviates from the usual Rayleigh distribution^ [43, 
72]. However, only few reports have mentioned this PDF asymmetry in laboratory 
experiments [74]. 

3.2. Intermittency in wave turbulence [35]. One of the most striking feature of 
turbulence is the occurrence of bursts of intense motion within more quiescent fluid 
flow. This generates an intermittent behavior [18, 57]. Indeed, a signal is called 
intermittent as soon as it displays more and more intense events when it is probed 
on a more and more short duration, leading thus to a strongly non Gaussian statis- 
tics at small time scales. The origin of non-Gaussian statistics in 3D hydrodynamic 
turbulence has been ascribed to the formation of coherent structures (strong vor- 
tices) since the early work of Batchelor and Townsend [18]. However, the physical 
mechanism of intermittency is still an open question that motivates a lot of studies 
in 3D turbulence [2, 19]. Intermittency has also been observed in a lot of prob- 
lems involving transport by a turbulent flow [passive scalar. Burgers turbulence of 
strongly compressible fluid] [39] , in magnetohydrodynamic turbulence in geophysics 
[26] or in the solar wind [1]. The possible observation of intermittency in wave 
turbulence that strongly differs from high Reynolds number hydrodynamic turbu- 
lence is of primary interest. It could indeed motivate explanations of intermittency 
different than the ones considering the dynamics of the Navier-Stokes equation. 

We have reported the first observation of intermittency in wave turbulence [35]. 
By measuring the temporal fluctuations of the surface wave amplitude, 'r]{t), at 
a given location (cf. Fig. 1), one can compute the local slope increments of the 
gravity surface waves at a certain time scale r: Arj{T) = rj{t + t) — 2rj{t) + rj{t — r). 
As shown in Fig. 6, starting from a roughly Gaussian statistical distribution of 
these increments at large scale t, the PDFs undergoes a continuous deformation 
towards strongly non-Gaussian shape at smaller scales. This shape deformation of 
the distributions across the scales is a direct signature of intermittency. Figure 7 
shows the structure functions of these increments, defined as 5p(r) = {\A7'i{t)\p) , 
which are found to be power laws of the time scale r on more than one decade, such 
as 5p(t) ^ where is an increasing function of the order p. The evolution of 
the exponents with p is found to be a nonlinear function of p as shown in Fig. 8. 
The nonlinearity of is another direct signature of intermittency. The slope at the 
origin of the curve ^p vs. p (the non-intermittent part) is close to ip/2, this value 
being easily derived by dimensional analysis'^. However, there exists any model that 
can describe the intermittent nonlinear part of at large p. 



Gaussian process with negative and positive values has a Rayleigh distribution for its crest- 
to-trough amplitude. 

^One has indeed 5p(t) ~ £P/6gP/V3p/2. 
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Figure 6. Probability density functions of the local slopes of wave 
amplitudes [■q{t + t) — 2ri{t) + r]{t — T)]/aT for different time lags 
6 < T < 100 ms (from top to bottom). Gaussian fit with zero 
mean and unit standard deviation (dashed line). Correlation time: 
Tc ~ 63 ms. Each curve has been shifted for clarity. 




4 10 20 30 40 50 

X (ms) 



Figure 7. Structure functions of order p of the local slopes of 
wave amplitudes, Spir), as a function of time lag t, for 1 < p < 6. 
(— ): Power-law fits, Sp{T) ^ t^p , where the slopes depend on 
the order p (see Fig. 8). Curves have been shifted for clarity. Inset: 
Flatness S4/S2 as a function of r fitted by a power law of slope 
—0.88 (— ). Correlation time: Tc — 63 ms. 



Recently, it has been proposed that theoretical corrections should be taken into 
account in weak turbulence to describe a possible intermittency that may be con- 
nected to singularities or coherent structures [20, 23] such as wave breaking [95] or 
whitecaps [23] on the fluid surface. However, intermittency in wave turbulence is 
often related to non-Gaussian statistics of low wave number Fourier amplitudes [20] , 
and thus it is not obviously related to small-scale intermittency of usual hydrody- 
namic turbulence. Here, we have shown the intermittent feature of the local slope 
of turbulent gravity waves in a same way as intermittency in usual turbulence. A 
challenge is to understand the physical origin of intermittency in wave turbulence. 
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Figure 8. Exponents of the structure functions as a function 
of the order p. is computed from the slopes of Fig. 7, and is 

fitted by ( ) = cip — with ci = 1.65 and C2 = 0.2. Solid 

line: dimensional analysis = 3p/2. 



3.3. Fluctuations of the energy flux [36, 37]. The key parameter of wave tur- 
bulence is the energy flux injected in the system. It is assumed to be conserved 
during the cascade across the scales. For the first time in wave turbulence, we have 
measured the instantaneous injected power in the fluid I{t) = F xV , where F{t) is 
the force applied by the wave maker on the fluid, and Vit) the wave maker velocity 
(see Fig. 1). The typical temporal evolution of this energy flux is shown in the 
inset of Fig. 9. It displays fluctuations much stronger than its mean value (/) (see 
Fig. 9). These fluctuations of energy flux have not been taken into account in the 
theoretical models that have been developed up to date. To wit, the knowledge of 
their evolution during the whole cascade would be necessary. Here, one studies the 
statistical properties of the energy flux fluctuations measured on the wave maker 
(integral scale). Numerous events of negative energy flux (/(t) < 0) are observed 
on Fig. 9 with a fairly high probability. These are events for which the wave field 
gives back energy to the wave maker. Understanding this phenomenon is necessary 
to develop marine renewable energy such as ocean wave energy. It is thus of first 
interest to model this empirical statistical distribution of the energy flux. To wit, 
the wave maker is described by a Langevin-type model with a pink noise (Ornstein- 
Uhlenbeck process) that mimics the experimental narrow- band random forcing [36] . 
Two linear coupled Langevin-type equation for V and F are then obtained; these 
two variables being Gaussian (due to the forcing) as observed experimentally (see 
Fig. 10). The computation of the distribution of the product of these two corre- 
lated Gaussian distributions thus gives the analytical expression of the energy flux 
distribution with no adjustable parameter [4, 36, 37]. Figure 9 shows a very good 
agreement between our model and the experiment. The asymmetry of the probabil- 
ity distribution is driven by the mean energy flux (/) , that is equal, in a stationary 
state, to the mean dissipated power. This model is not restricted to wave turbu- 
lence, but is very generic since it also describes the energy flux distribution in other 
dissipative out-of-equilibrium systems such as turbulent convection, granular gases, 
bouncing ball with a random vibration [4], and electronic circuit with a random 
voltage (see below and Ref. [37]). 
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Figure 9. Inset: Time recording of injected power, I{t), driven the 
surface waves. (J) = 2 niW. Main: Probability density function 
of injected power showing two exponential tails. Its asymmetry 
is related to its mean value (/) (vertical solid line). Note that 
fluctuations are much more larger the mean value. Blue dashed 
line: prediction of our model with any adjustable parameter. 
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Figure 10. Time recordings of the velocity V{t) of the wave maker 
(left) and the force applied F{t) to the wave maker by the vibration 
exciter (right). Both probability distributions are quasi Gaussian 
(see dashed lines) with zero mean value, ay and ap are the stan- 
dard deviations. 



We have also emphasized the bias that can result from the system inertia when 
one tries a direct measurement of the fluctuations of injected power. When the wave 
maker inertia is not negligible^, the measurement of the instantaneous force has to 
be corrected by the inertia before performing the product with the velocity signal 
in order to deduce the power. Otherwise, it may lead to wrong estimations of the 
fluctuations of injected power, the mean value being not obviously affected. There 
exists only a few previous direct measurements of injected power in 3D turbulent 
flows, and this type of inertial bias has never been taken into account [61, 76, 89]. 



For our experiments, this is the case with water but not with mercury. 
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Frequency (Hz) 

Figure 1 1 . Spectra of the surface wave amplitude divided by the 
variance ay of the wave maker velocity for different forcing ampli- 
tudes: av = 2.1,2.6,3.5 et 4.1 cm/s. The dashed line has slope 
—5.5, whereas the solid line has slope —17/6. Inset : The mean 
injected power (/) is found proportional to ay. 

The mean injected power (/) is found to be proportional to the fluid density^'^ 
p and to ay, the rms value squared of the wave- maker velocity fluctuations. The 
rms value of the injected power fluctuations also scales as aj ^ (/) ~ pcry. The 
energy transfer occurs via 4- wave interactions in the gravity regime and via 3-wave 
interactions in the capillary regime. Thus, if these two regimes are independent as 
in Eqs. (5) and (6), the dependence of the spectrum amplitude Sri on the mean 
energy flux (/) should display different scaling laws. Our measurements show that 
a single scaling law, Sri ~ (/), leads to the collapse on a single master curve of our 
experimental spectra for different values of the forcing (see Fig. 11). This strong 
disagreement with the weak turbulence theory [cf. Eqs. (5) & (6)] thus shows that 
interactions between these regimes are very important although they are not taken 
into account theoretically. 

The experimental distribution of the injected power averaged on a time lag r, 
P{It), has also been studied [36] in the framework of the Fluctuation Theorem 
of the out-of-equilibrium statistical physics (also called the Gallavotti-Cohen rela- 
tionship) [33, 44]. This theorem states that - e-/3(^)^-/(-f> where P[x] 
is the probability to have the value x and /3(t) is a constant related to the en- 
ergy of the system. The theorem hypotheses (such as the temporal reversibility) 
are generally not fulfilled for a real dissipative system. The apparent verification 
of this theorem in numerous previous experiments [21, 22, 41, 85, 86] is due to 
their small range of explored fluctuation amplitudes It/ {I) at long averaging time 
T [3, 5, 6, 45, 77, 91]. Our wave turbulence experiment allows us to reach large 
enough fluctuations of injected power, and we have shown that they do not verify 
the theorem [36], but are found in good agreement with the analytical calculation 
performed with a Langevin-type equation forced with a Gaussian white noise [40] . 

We have also studied experimentally the energy flux fluctuations in an electronic 
RC circuit driven with a stochastic voltage [37]. This out-of-equilibrium model 
system allows us to easily control the system parameters (notably, its damping rate 



^"This has been verified for only two different fluids (water and mercury) but with a density 
ratio close to 14. 
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Figure 12. Fluctuations of injected power I{t) in an electronic 
RC circuit driven with a stochastic voltage for two different values 
of the damping rate {RG)~^ = 2000 and 200 Hz (resp. blue and 
red curves). Comparison between experiment (— ) and theoretical 
prediction of the Langevin-type model ( ). As in wave turbu- 
lence, the more dissipative the system is, the more asymmetric the 
PDF of I{t) is. 

l/(i?C)). We show that the large fluctuations reached do not verify the Fluctuation 
Theorem even for long averaging time r. As shown in Fig. 12, the probability 
distribution of the instantaneous energy flux is qualitatively similar to the one of 
wave turbulence of Fig. 9, and is indeed well described by our above model (see 
dashed lines in Fig. 12). This electronic circuit is thus one of the simplest system 
to understand certain properties of the energy flux fluctuations shared by other 
dissipative out-of-equilibrium systems, such as in granular gases [41], convection 
[85, 86], or wave turbulence [36]. 

3.4. Wave turbulence in low-gravity environment [38]. Experiments of wave 
turbulence in low gravity are relevant to study pure capillary waves on a large 
range of wavelengths without interfering with the gravity waves. During ground 
experiments, the upper limit is the capillary length, and its lower limit is related to 
viscous dissipation. Our experiment in weightlessness has been carried out during 
CNES or ESA campaign of parabolic flights on board of an specially modified 
Airbus A300. A campaign consists of 3 fiights, each flight consists of a series of 30 
parabolic trajectories that result in low-gravity periods, each of 22 s. Such a low- 
gravity environment (about ±5 x 10~^g) allows us to also study the behavior of a 
spherical fluid layer on which waves can be propagates with no boundary condition 
(see the experimental setup in Fig. 13). 

In such low-gravity conditions, we have observed the capillary wave turbulence on 
the surface of a fluid layer that covers all the internal surface of a spherical container 
which is submitted to random forcing [38]. Figure 14 displays the capillary wave 
spectrum over a large range of wavelengths, notably in wavelengths where gravity 
waves were present during ground experiments. This spectrum over two decades in 
frequency (corresponding to wavelength from mm to few cm) is found in roughly 
good agreement with the frequency-exponent prediction of weak turbulence theory 
of the capillary regime in —17/6 ~ —2.8 [see Eq. (6)]. As shown in Fig. 14, this 
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Figure 13. Experimental setup for the wave turbulence in low 
gravity. In low-gravity periods, the fluid covers all the internal 
surface of the tank which is submitted to vibrations. 




Figure 14. Power spectra of surface wave amplitudes in low gravity 
for two different forcings: random forcing from to 6 Hz (lower), 
or sinusoidal forcing at 3 Hz (upper). Dashed lines have slopes of 
-3.1 (lower) et -3.2 (upper). Inset : Same with gravity. Two power 
laws appear corresponding to the gravity (slope -5) and capillary 
(slope -3) regimes. 



exponent is also found to be independent of the large-scale forcing parameters (low- 
frequency sinusoidal or random forcing) as expected within the inertial range. 

3.5. Wave turbulence on the surface of a ferrofluid [11]. A fcrrofluid is a 
colloidal suspension of nanometric ferromagnetic particles diluted in a liquid dis- 
playing thus both liquid and strongly magnetic feature. In contrast with usual 
liquids, the dispersion relation of surface waves on a ferrofluid is nonmonotonic and 
displays a minimum (depending on the amplitude of the applied magnetic field B 
[16, 17, 32, 65, 69]) for which an instability occurs. Above a critical field Be, a static 
hexagonal pattern of peaks occurs on the ferrofluid surface, the so-called Rosensweig 
instability [25]. The amplitude of the magnetic field, B, can thus easily tune the 
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Figure 15. Experimental setup for the study of wave turbulence 
on the surface of a ferrofluid in a normal magnetic field generated 
by two coaxial coils. Surface waves on the ferrofluid are generated 
by the horizontal motion of a plunging wave maker driven by an vi- 
bration exciter. The amplitude of surface waves at a given location 
is measured by a capacitive wire gauge. 



dispersion relation of the surface wave which reads [13, 81] 

where fiQ ~ An x 10"'' H/m is the magnetic permeability of the vacuum, and f[x] 
a known function of the magnetic susceptibility of the ferrofluid^^. One can thus 
easily tune, with just one single control parameter i?, the dispersion relation of 
surface wave from a dispersive one (fc or term) to a nondispersive one (fc^ term). 
Besides, some theoretical questions arc open notably about the possible existence of 
power-law solution for the energy wave spectrum for two-dimensional nondispersive 
systems [23, 63, 71, 102]. Wave turbulence on the surface of a magnetic fluid should 
thus display striking properties depending on whether the system is below, close or 
above the onset of the Roscnswcig instability. 

We have studied the wave turbulence on the surface of a ferrofluid submitted to a 
normal magnetic field (see Fig. 15) [11]. We have shown that magnetic surface waves 
arise only above a critical field, but below the onset Be of the Rosensweig instability. 
As shown in Fig. 16, the power spectrum of the magnetic wave amplitudes displays 
a single power law over the whole range of accessible frequencies instead of two 
regimes (gravity and capillary) observed with no applied magnetic field. On can 
compute by dimensional analysis the power spectrum of magnetic wave turbulence 
as 

5,"'^^.^) ^ f — ) ' u-^ , (8) 
\PP-oJ 

This predicted frequency exponent uj~^ is too close to the one predicted for the 
capillary regime in to be distinguished experimentally with our experimen- 

tal accuracy. This could explain that only one single slope is observed in Fig. 16 
which thus corresponds to a "magneto-capillary" wave turbulence regime. Unlike 
the above dispersive systems [see Eq. (5) & (6) and §3.1], this —3 exponent does 



^^Indeed, x(H) depends on the applied magnetic field, H, through Langevin's classical law, 
and thus on the magnetic induction B, since B = fio{l + x)^- The magnetic induction, B, will 
be called hereafter the magnetic field for simplicity. 
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Figure 16. Power spectrum of wave amplitude, ij{t), on the fer- 
rofluid surface for two applied magnetic field. Inset: B = : Grav- 
ity and capillary wave turbulence regimes. Dashed lines have slopes 
-4.6 and -2.9. Vertical dashed line: crossover frequency fgc Hz. 
Main: B = O.QBc : Magnetoeapillary wave turbulence. Dashed 
line has slope -3.3. Crossover is decreased towards fgc ~ 5 Hz. 



not depend on the energy flux exponent a. However, the exponent a is still not 
deduced by dimensional analysis. By measuring the scaling law between the spec- 
trum amplitude 5^"^ and the magnetic field B, one then shows that a = 1/3, i.e. 
that the magnetic wave turbulence involves a 4-wave interaction process. One has 
thus obtain the whole expression of the spectrum of magnetic wave turbulence. 




Figure 17. Crossover frequencies between gravity, magnetic and 
capillary wave turbulence regimes as a function of the applied mag- 
netic field. Different forcing parameters: (x) 1 to 4 Hz, (o) 1 to 5 
Hz, and (+ or o) 1 to 6 Hz. Theoretical curves fgc, fgm and fmc 
can be computed analytically from Eq. (7) [11]. The triple point 
reads ft = ^{^f^ - 10.8 Hz and = / f[x{Ht)], that 

is Bt/Bc = 0.65 for our ferrofluid, where Be is the onset of the 
Rosensweig instability. 

The existence of the regimes of gravity, magnetic and capillary wave turbulence 
is reported in the phase space parameters as well as a triple point of coexistence of 
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these three regimes (see Fig. 17) [11]. The theoretical crossover frequency curves, 
as well as the triple point, can be computed from the dispersion relation of Eq. (7), 
and arc found in good agreement with the experimental data with no adjustable 
parameter (see solid lines in Fig. 17). Finally, the statistical distribution of the 
wave amplitudes and their change with the magnetic field has been also studied; 
our measurements confirm the nonlinear nature of the wave interaction in a wave 
turbulence regime, and also show the onset of existence of the magnetic waves [11]. 

4. Conclusions and open problems. I have first presented the context and the 
state of the art in wave turbulence. I have then focused on the wave turbulence 
on the surface of a fiuid, and I have discussed the main experimental results of our 
group on this subject. We have measured, for the first time, the key governing 
parameter of wave turbulence: the energy fiux that drives the waves and cascades 
to small scales through nonlinear interactions. Fluctuations much larger than the 
mean value have been observed as well as instantaneous negative energy fiux events 
that occurs with a fairly large probability. Taking into account these energy flux 
fluctuations in theoretical models of cascades remains an open problem. The scaling 
of the power spectra of capillary wave turbulence is found in agreement with weak 
turbulence theory both in laboratory experiments or in low-gravity environment 
(to study pure capillary waves). A disagreement with weak turbulence theory has 
been observed for the gravity wave spectrum. It has been also found that gravity 
waves are intermittent, i.e. the statistics of velocity increments of the interface 
are strongly non Gaussian at small scales. This flrst observation of intermittency in 
wave turbulence could be related to the large fluctuations of the energy flux. Finally, 
magnetic wave turbulence has been observed on the surface of a magnetic fluid 
(ferrofluid) as well as a triple point of coexistence of gravity, capillary and magnetic 
wave turbulence regimes. These features are understood using dimensional analysis, 
but no weak turbulence prediction exists for this magnetic case. 

As emphasized in this review, there has been a renewal of interest in wave tur- 
bulence since last years. However, open problems are still numerous. It is crucial to 
consider in a precise way the range of validity of theoretical predictions and to study 
experimental effects associated to breakdown of existing theories. Finite size effects 
in experiments lead to a quantization of wave number that can strongly affect en- 
ergy transfers. Understanding the observed dependence of the frequency-exponent 
of the gravity wave spectrum with the forcing parameters is thus of primary inter- 
est. Besides, most in situ or laboratory measurements on wave turbulence are local 
in space whereas theoretical predictions often concern Fourier space, e.g. power 
spectra or probability density functions of the fluctuations of the Fourier amplitude 
of a wave component at a given wave number. Thus, an important challenge is to 
develop a direct probe in Fourier space. Some recent studies measure the ampli- 
tudes of gentle interacting waves on a spatial zone [67, 87, 93, 94]. However, these 
optical methods cannot measure amplitudes of steep nonlinear waves as involved 
in wave turbulence. Nevertheless, one need to access to such steep wave amplitude 
in the Fourier space in order to get a better understanding of the elementary dy- 
namical processes involved in the energy cascade. An other challenge is to observe 
a possible inverse cascade in wave turbulence that has been predicted theoretically 
[97, 99] and numerically [58]. Finally, a significant issue is to understand the physi- 
cal mechanisms of intermittency in wave turbulence that has been recently observed 
[35]. 
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